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. a first dynamical calculation of the widths of the 0+ and 8* pentaquarks for finite Nc in a 

' chiral soliton model. 

(N 

Q , PACS numbers: 12.39. Dc, 13.30.Eg, 13.75.Jz 

U ' Keywords: chiral solitons, pentaquarks, kaon-nucleon scattering, resonance phase shift 

Q 

; I. INTRODUCTION 

• Chiral soliton model predictions for the mass of the hghtest exotic pentaquark, the 0"*" with 

l/^ ! zero isospin and unit strangeness, have been around for about two decades P, 0, 0, H, 1^- 

' Nevertheless, the study of such pentaquarks as potential baryon resonances became popular only a 

■ short while ago when experiments indicated their existence^. These experiments were stimulated 

ly^ , by a chiral soliton model estimate that claimed that the width of the pentaquark decaying into 

I a kaon and a nucleon would be much smaller than those typical for hadronic decays of baryon 
resonances 0, U, however, see also 0]. Such narrow resonances could have escaped detection 



in earlier analyses. 

i-^ ■ A. Hadronic Decays in Soliton Models 



^ , Estimates for the G"^ width are based on identifying the interpolating kaon field with the 

$H ' (divergence) of the axial current, which in turn is computed solely from the classical soliton. This 

is essentially a generalization to flavor SU (3) of the early calculation of the A — > vrA^ decay width 
in the pioneering paper of ref. in the two flavor Skyrme model. Since then there have been 
quite some efforts to understand that process in chiral soliton models by going beyond the sirnple 
identification of the pseudoscalar meson fields with the axial current [lil . Efi. .17,, 13,, 1^ . 0, 2^, 



l23l |: see also chapter IV in ref. |2J]. It must be stressed that the description of hadronic decays 
in soliton models is a severe problem and has become known as the Yukawa problem in soliton 
models. The reason is that the interaction Hamiltonian that mediates such hadronic transitions 
must be linear in the field operator that corresponds to the pseudoscalar meson in the final state. 
On the other hand the naive definition of the soliton as a stationary point of the action prohibits 
such a linear term. Hence in these models resonance widths must be extracted from appropriate 
S'-matrix elements of meson baryon scattering. Most of the studies on meson baryon scattering 
are restricted to the adiabatic approximation [2^, E^ . 2^ . 2^ which represents the next -to-leading 



order in the 1/Nc expansion where Nc is the number of colors. Unfortunately in the large- A'^c' 



^ For recent reviews on the experimental situation we refer to refs. [aH- The numerous theory papers that have ap- 
peared since then may be downloaded from jhttp: //www. r cup . osaka-u. ac . jp/~hyodo/research/Thetapub.html| 
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limit the A is degenerate with the nucleon and its width is of even lower order in that expansion. 
Therefore its width cannot be computed in adiabatic approximation and the above cited studies 
mainly motivate Born contributions to the ^-matrix element of pion nucleon scattering in the P33 
channel at 0{\/Nc)- For the pentaquarks the situation is very different. As we will see, the 0"*"- 
nucleon mass difference is non-zero even as Nc 00 and hence its hadronic decay width should 
be at the same order as the adiabatic approximation for kaon-nucleon scattering. In that respect 
the hadronic decays of pentaquarks are a perfect playground as any treatment that induces Born 
contributions can be tested against the adiabatic approximation. 



B. Quantizing the Soliton in Flavor SU(3) 

When treating chiral solitons in flavor SU{3) we face the problem how to properly generate 
baryon states from the classical soliton. Essentially there are two methods at our disposal, the 
bound state approach (BSA) [2^ and the rigid rotator approach (RRA) js^]; see e.g. ref. [sj] for a 
review and further references. The latter is essentially the generalization of the collective coordinate 



quantization adopted earlier for the two flavor model [ij]. In this approach collective coordinates 
that describe the spin flavor orientation of the soliton are introduced and quantized canonically. As 
this corresponds to quantizing a rigid top in SU (3) it is natural that the corresponding eigenstates 
are characterized according to SU (3) multiplets, most prominently the octet (8) and decuplet 
(10) for baryons with spin J = ^ and J = |, respectively. However, also more exotic SU{3) 
representations like the anti-decuplet (10) emerge. In addition to states with quantum numbers of 
octet baryons this representation also contains states with quantum numbers that cannot be built 
as three-quark composites. These states contain (at least) one additional quark-antiquark pair 
which lead to the notion of exotic pentaquarks. Most prominently there have been experimental 
claims 05 0) 0| for a pentaquark, the strangeness S = +1 isosinglet that sits at the top of the 
anti-decuplet. 

In the BSA kaon fluctuations i]a{x,t), with a = 4, . . . , 7, are introduced and quantized as 
harmonic oscillations. It corresponds to the adiabatic approximation and should be exact as Nq 
00. Common to all soliton models is the appearance of a bound state with strangeness S = —1 
which serves to construct states that correspond to the ordinary hyperons. A, S and S*. On the 
other hand neither a bound state nor a clear resonance has been observed in the 5 = +1 channel. 
This non-observation has often been used to argue that the pentaquark prediction would be a mere 
artifact of the RRA Here we will show that this statement is completely based on a 

misinterpretation of the RRA. For a sensible and ultimate comparison of the BSA and the RRA 
it is necessary to also include small amplitude fluctuations in the RRA. We call that approach 
the rotation-vibration approach (RVA). However, the fluctuations in the RVA must be constrained 
to the subspace that is orthogonal to the rigid rotation. We therefore label these fluctuation by 
f]^{x, t). Then the test mentioned at the end of subsection A is to first compute the (S-matrix from 
the constrained fluctuations, r]^{x,t) and supplement that result with the contributions from the 
(iterated) O"*" exchange. That yields the wave-functions rjaix,t) that are constrained to the same 
subspace as fj^{x , t). The rjaix ,t) are the full solution to the RVA. In the limit Nc co rja{x,t) 
must yield the same (S-matrix as in the BSA, i.e. when extracted from r]a{x,t). In performing 
this test, our studies serve three major purposes: First, we will establish the equivalence of BSA 
and RVA in the large-iVc limit. Second, we will compare ?7q, and r]^ to extract a distinct S = 1 
resonance in kaon nucleon scattering that survives as Nq — > 00. Third, we will provide the first 
dynamical calculation of the width in a chiral soliton model for finite Nc. 
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C. General Remarks 

Let us round off this introduction section with a few general remarks on the comparison between 
the BSA, the RRA and the RVA. Generally we may always introduce collective coordinates to 
investigate specific excitations. This is irrespective of whether the corresponding excitation energies 
are suppressed in the large-A'^c' limit or not. That is, in contrast to the claims of refs. [s^] no 
separation of scales is required to validate the collective coordinate quantization. Eventually the 
coupling terms and constraints ensure correct results. Even more, we will observe in section V 
that collective coordinates may also be introduced for mx 7^ w-^r- This result generalizes to the 
statement that collective coordinates are not necessarily linked to the zero-modes of a system 
and any distinction between dynamical and static zero-modes js^] is obsolete. Already some time 
ago the concept of describing baryon resonances as collective excitations has been successfully 
applied to quadrupole [s^ and monopole excitations, both in two and three [s^l flavor models. 
In particular collective monopole excitations have been studied in the context of the pentaquark 
excitation 3^, 3^ . 

Sometimes the correct description of the A(1405) as an 5'-wave bound state is considered as 
an argument in favor of the BSA vs. the RRA. Again it must be stressed that the full RVA also 
contains kaon fluctuations and thus also describes the formation of an S-wave bound state. As 
shown in ref. this properly accounts for the A(1405) in the RVA. 



D. Outline 

This paper is organized as follows. In section II we will deflne the model and describe the fleld 
parameterization in the bound state and rigid rotator approaches in more detail. In sections III 
and IV we extensively compare these approaches and show how they yield identical spectra in the 
limit Nc — > 00, even for non-zero symmetry breaking, when the fluctuating modes are limited 
to the same subspace. In section III we also show that the constraints for r]^ serve to separate 
resonance and background pieces of the phase shift. Section V represents a main piece of our 
paper as we combine rotational and fluctuation degrees of freedom to construct the solution for 
the RVA up to quadratic order in the fluctuations. We employ this solution to compute the width 
of the G"*" pentaquark for rriK = rn-,^ in section VI. Another major result is that in the flavor 
symmetric case only a single collective coordinate operator mediates the transition — > KN at 
the leading order of the large-A^c expansion. This is in direct contradiction to approaches that 
construct transition interactions with (at least) two SU{2>) operators and adjust the coefficients to 
yield a narrow 6+ 0, fill . I41I . |43 | . In section VII we extend our treatment to the physical case 
rriK 7^ w-TT- This brings a second SU{3) operator into the game. In section VIII we collect our 
numerical results and summarize in section IX. We relegate technical details on the derivation of 
the equation of motion for the constrained fluctuations to an appendix. 



II. THE MODEL 

For simplicity we consider the Skyrme model [i^ as a particular example for chiral soliton 
models. However, we stress that our qualitative results do indeed generalize to all chiral soliton 
models because these results solely originate from the treatment of the model degrees of freedom. 

A preliminary remark on notation: In what follows we adopt the convention that repeated 
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indices are summed over in the range 



a,b,c, . . . 
a,/3,7, . . 
i, j, k, . . . 



1,...,8 
4,. ..,7 

1,2,3 



(kaonic) 
(pionic) . 



(2.1) 



Chiral soUton models are functionals of the chiral field, U, the non-linear realization of the 
pseudoscalar mesons, (pa 



U{x ,t) = exp 



-(/>a(£,t)Aa 



(2.2) 



with Aa being the Gell-Mann matrices of SU{3). For a convenient presentation of the model we 
split the action into three pieces 



r = TsK + ^wz + ^SB ■ 
The first term represents the Skyrme model action 

1 



SK 



+ 



32e^ 



(2.3) 



(2.4) 



Here /jr = 93MeV is the pion decay constant and e is the dimensionless Skyrme parameter. In 
principle this is a free model parameter. The two-flavor version of the Skyrme model suggests to 
put e = 4.25 from reproducing the A-nucleon mass difference. The anomaly structure of QCD is 
incorporated via the Wess-Zumino action [3 | 



wz 



iNc 
2407r2 



(2.5) 



with = U^df^U. Note that T\yz vanishes in the two-flavor version of the model. Furthermore 
Nc is the number of colors, the hidden expansion parameter of QCD 4^, 4^ . 
The flavor symmetry breaking terms are contained in T^b 



d'^xii m(u + U^ -2 



with 



M 



(2.6) 



2m\- 



Unless otherwise stated we will adopt the empirical data for the masses of the pseudoscalar mesons: 
ttItt = 138MeV and mx = 495MeV. We do not include terms that distinguish between pion and 
kaon decay constants even though they differ by about 20% empirically. This is a mere matter 
of convenience rather than neg ligence. This omission leads to a considerable underestimation of 
symmetry breaking effects |47| which approximately can be accounted for by rescaling the kaon 
mass niK rnKfR/fn- 

Although physically sensible results are obtained by setting Nq = 3, it is very illuminating 
to organize the calculation in powers of 1/Nc and consider Nc undetermined. In the context of 
this Nc counting the mass parameters are 0{N^) while f-^ = 0{yJNc) and e = 0{1/ \/Nc) to 



ensure that the (perturbative) n-point functions scale as A'^^ |46l |. In order to study the Nc 
dependence in the three flavor version of the soliton model we include the Nc dependence in the 
choice of the model parameters via 



U = 93MeV 



Nc 



and 



4.25 



Nc 



(2.7) 
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The action, eq. ()2.3() allows for a topologically non-trivial classical solution, the famous hedgehog 
soliton 



Uo{x 



exp 



iX ■ xF{r) 



\x 



(2. 



which is embedded in the isospin subspace of flavor 5J7(3) parameterized by the Gell-Mann ma- 
trices Aj, i = 1,2,3. Other embeddings are possible but energetically disfavored due to flavor 
symmetry breaking. The chiral angle, F{r) is obtained from integrating the classical equation of 
motion extracted from eqs. (|2.4j) and (|2.6|) subject to the boundary condition^ -^(0) — F{oo) = vr 
ensuring unit winding (baryon) number. The soliton can be constructed as a function of the 
dimensionless variable efT^r and is thus not subject to Nc scaling. 

The excitations of this soliton configuration are quantized to generate physical baryon states. 
The central task in extending this picture to three flavors is the incorporation of strange degrees 
of freedom, that is kaons. In a first step we introduce fluctuations^ r]a{x,t) with a = 1, ... ,8, for 



the pseudoscalar fields 2^, 4^, 



U{x,t) = A/[/o(^)exp 



— XaVaiXjt) 
Jtt 



(2.9) 



Expanding the action in powers of these fluctuations actually is an expansion in rja/ fn and thus a 
systematic series in I/^/Nq- The term quadratic in rja is of special interest. It describes scattering 
of mesons off a potential that is generated by the classical soliton, eq. (|2.8|) . In particular an 
anti-kaon is bound in the P-wave channel. Combined with the soliton it may be interpreted as the 
A hyperon. However, the system consisting of the soliton and this bound anti-kaon has still to be 
projected onto states with good spin and isospin. This is accomplished by introducing collective 
coordinates for the spin-isospin orientation^, ^2(t) G SU{2) 



Uix,t)=A2it)^/Uo{£] 



exp 



-Xaria{x,t) 



VUoW)Alit) 



(2.10) 



and quantizing them canonically 2^, 4^ . In this treatment the kaon fluctuations a = a = A, ... ,7 
decouple in the adiabatic equations of motion. This leads to the so-called bound state approach 
(BSA) that we will review in section III. 

The seemingly alternative approach is to allow the collective coordinates to span the full fla- 
vor space, ^3(i) G SU{3) and quantize A3 canonically to generate baryon states with non-zero 
strangeness while omitting the fluctuations. Baryon states are then classified according to SU{3) 
representations. This is the so-called rigid rotator approach (RRA) [sO]- However, for a sensi- 
ble investigation of the rotation-vibration coupling the collective coordinate approach must be 
completed to also contain meson fluctuations. Hence the ansatz that defines the RVA reads 



U{x,t)=As{t)y'Uoix) 



exp 



— Xar]aix,t) 
In 



y^Uo{x)As{t)K 



(2.11) 



Obviously eq. (|2.11|) comprises the ansatz, eq. (|2.10j) . In contrast to eq. (|2.10j) the ansatz (|2.11|) has 
additional collective rotations into strange directions. But these field components are also contained 
in the fiuctuations, rja. It is therefore necessary to constrain these fiuctuations to the subspace that 



^ The overall sign is a matter of convention. 

^ In what follows we will discuss different parameterizations of the chiral field, U. We abstain from introducing 

additional labels to distinguish among them. 
* Due to the hedgehog symmetry only one set of collective coordinates is required. 
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is orthogonal to the rigidly rotating hedgehog. These constraints modify the equations of motion 
for the fluctuations, such that rja 7^ r]a- This is not new: In principle we encounter the same 
problem already in SU (2) where the ansdtze eqs. ()2.1U() and 1)2.11(1 are identical. The difference 
to the SU (3) case is that the redundant modes are automatically orthogonal to the rigid rotations 
when treating the SU{2) model in adiabatic approximation. This changes, however, as soon as 
non-adiabatic contributions or external fields e.g. photons are considered. We will discuss 
and derive the modifications that occur in the SU{3) case in section III and a separate appendix, 
respectively. 

It is obvious that the ansdtze (j2.1(H) and (j2.11j) span the same Fock space in the baryon number 
one sector. An essential difference, however, is that an expansion to a given order in rja using 
eq. ()2.1()|) contains the kaon fields only up to a fixed power in the I/Nq expansion while the 
ansatz (|2.11() treats the collective rotations to all orders. Obviously the two treatments must yield 
identical results in the limit Nq — > 00. This provides a thorough check on our calculations. We 
also stress that the RVA formulation induces terms in the action that are linear in rja (in contrast 
the BSA does not give terms linear in r]a)- This gives rise to Yukawa couplings from which we will 
evaluate widths of rotational excitations. 

There have been earlier studies of the non-adiabatic rotation-vibration coupling in the Skyrme 
model. In the SU{2) version the contribution linear in the time derivative of A2 has been shown to 



properly describe the Tomozawa- Weinberg relations for S'-wave pion nucleon scattering |5ll | . Also 
in the two flavor model, the terms quadratic in the time derivative of the collective coordinates 
have been employed to incorporate thresholds and A resonance exchange effects in pion-nucleon 
scattering |2nl. Ell. E3. . In the three-flavor BSA the rotation-vibration coupling gives rise to 
the hyperfine mass splitting [i^ that we will discuss in the following section. Also in SU (3) meson 
baryon scattering has been investigated in the RVA Unfortunately, the strangeness, S = +1 
channel, which contains the pentaquark, was not considered in that outstanding paper. 

III. BOUND STATE APPROACH (BSA) 

In this section we will discuss the construction of P-wave baryons with strangeness = ±1 
in the framework of the BSA starting from the parameterization, eq. (|2.9|1 . Baryons with zero 
strangeness are treated as in two-flavor chiral soliton models (l^ . 

A. Conventional BSA in the P wave channel 

As in refs. [2^, li^ we expand the action to quadratic order in the kaonic fluctuations rja in 
the background of the soliton, eq. (|2.8|) . These wave-functions describe kaon-nucleon scattering 
in the intrinsic frame. In this frame the modes decouple with respect to grand spin Q = i ziz ^, the 
eigenvalue of the sum of kaon orbital angular momentum (i) and isospin. This defines the intrinsic 
T-matrix, Tcp . The T-matrix elements in the laboratory frame T{£I2J) are linear combinations 
of the Tg^£ 2a, |53|. For the G"*" channel (I = 0, J = ^) this linear relation is very simple 



r(POl) = Ti^ (3.1) 



2 

i.e. only the intrinsic P-wave kaon with grand spin G = ^ contributes. Its wave-function is 
characterized by a single radial function, that parametrically depends on the frequency in a Fourier 
analysis 



'TV^') i^'t) = r d^^'^'vUr)£-rx{c^)- (3.2) 
V% + ^V7j p J -00 
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Here xi'^) is a two-component iso-spinor that contains the amphtudes of the kaon mode. Upon 
quantization its components will eventually be elevated to creation- and annihilation operators. 
Of course, an analogous decomposition exists for rja that defines the radial function rji^{r). In what 
follows we will omit the subscript lo. The amplitude x drops out from the equation of motion that 
can be formulated for the radial functions, 



n(r) 

Mx(r) 



[ojMK{r) - A(r)] r/(r) 
A(r) 



UJ 



MK{r) 



n(r) 



(3.3) 



Up to a factor i the radial function n(r) parameterizes the momentum conjugate to the P-wave 



fluctuation j^. The radial functions in eq. (|3.3|) can be taken from the literature [IS*. 
terms of the chiral angle F = F{r), 



m 



A(r) 



hi + 



Nc wSin^F 



m: 



I with 



Mxir) 



1 + 



1 



+ 2 



sin^F 



dr"^ r dr ° ' 



(3.4) 



where F' = . The effective potential Vcsir) is quite an involved function that also contains 
spatial derivative operators. Since it does not play any special role for our studies we do not repeat 
it here. The radial function A(r), proportional to the baryon density, -Bo('')) originates from the 
WZ-term, eq. (|2.5() and thus contains the explicit factor Nc- The radial function Mk{t) is the 
metric function for kaon fluctuations. This becomes obvious when eliminating the momentum Ii[r) 
in favor of a second order differential equation for 7/(r), 



h^ T]{r) + LO [2A(r) - wM/^(r)] r/(r) 



0. 



(3.5) 



as Mfcir) multiplies the term of highest power in the frequency, to. Eq. (|3.5() could have been 
easily obtained without reference to the conjugate momentum. However, we have introduced it 
here because it becomes essential when implementing the constraints that characterize rj. The 
appearance of factors Nc/f^ and /^e^ in the radial functions shows that r/(r) does not scale with 
Nc- Hence terms in the action that are of higher than quadratic order in rja are suppressed due 
to the additional factor l/fn- Thus solving eq. (|3.5|) provides the exact T-matrix in the limit 
Nc oo. 

The equation of motion l|3.5j) is not invariant under particle conjugation lo ^ —uj, yielding 
different results for kaons {lo > 0) and anti-kaons (u < 0). This difference obviously originates 
from the Wess-Zumino term. It is well known that equation (|3.5)) has a bound state solution at 
io = —iOA, that gives the mass difference between the A-hyperon and the nucleon in the large-A'^c' 
limit. As this energy eigenvalue is negative it corresponds to a kaon, i.e. it carries strangeness 
S = —1. In the symmetric case {mx = rn-,^) the bound state energy vanishes, i.e. the bound 
state is nothing but the zero mode of SU{2>) flavor symmetry. The radial function of the properly 
normalized zero mode reads 



z{r) 



: sm- 



with Qk = fl j d^rMKir) sin'^^^ = 1.894GeV-i (^^^ . (3.6) 



It satisfies the important relation /iq z(r) = 0. The normalization constant, @k is the moment of 
inertia for flavor rotations into strangeness direction. In what follows we will employ the expressions 
"zero-mode" and "rotational mode" synonymously. 
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When symmetry breaking is switched on, the bound state evolves continuously from this zero 
mode. For uik = 495MeV we find = 196MeV. This bound state energy can be reliably 
estimated by sandwiching the equation of motion H3.5|) between zero mode wave-functions and 
solving for uj. We first define 



= J r'^drz'^{r)2\{r) 



Nc 
49^ 



0.396GeV 



and 



r = — — [mj^ - I r 



drz^ir) = 2.969GeV-^ (m|- - ml) 



Nc 



which are 0{N^) and 0{Nc), respectively. Then the equation for uj reads 



3T 
86^ 



+ UJOUJ . 



This quadratic equation for uj has the negative valued solution —uj\ with 



UJA 



3r 

26^ 



(3.7) 



(3.9) 



(3.10) 



For mK = iT^TT we have P = and thus also = from this estimate. In the realistic case we 
find r = 671MeV (^^) yielding UJA ~ 217MeV which is only about 10% larger than the energy 

eigenvalue of the exact solution to eq. (|3.5|) . Obviously the energy eigenvalue wa is a non-linear 
function of flavor symmetry breaking. In particular the observed relation 3r/20i^ > uJq indicates 
that the perturbative treatment of flavor symmetry breaking is bound to fail^. 

The quadratic equation (|3.9|) has a second solution: ujq. It corresponds to strangeness = +1 
fluctuations and determines the mass difference between the 0"*" pentaquark and the nucleon, 



1 



3r 

26^ 



+ Wo 



(3.11) 



Also this eigenmode evolves from the zero mode. However, the Wess-Zumino terms shifts ujq 
to the positive continuum. We therefore expect a resonance structure in the phase shift around 
ujQ = 396MeV and 617MeV for the cases mx = m-,^ and mx / 'tj-td respectively. These phase 
shifts are shown in figure ^ Apparently no clear resonance structure is visible; the phase shifts 
hardly reach 7r/2 at the estimated energies. The absence of such a resonance has previously lead 
to the conjecture and criticism that there would not exist a bound pentaquark in the large-A'^c 
limit [s^] and that it would be a mere artifact of the RRA. The absence of pentaquarks in the BSA 
does not seem to be a result limited to the simple Skyrme model. For example, in the BSA to the 
NJL chiral quark soliton model a pentaquark bound state has not been observed either [s^. 

We have obtained the above equations for wa and ujq by projecting the BSA equation (|3.5jl 
onto the rotational zero mode. Hence they are exact when the fluctuations are constrained to the 
rotational subspace. 



^ The realistic case is even worse, as F is underestimated by about 50% in our approximation with /k = f-w- 
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BOUND STATE APPROACH 



BOUND STATE APPROACH 
mK = 495MeV 



FIG. 1: The phase shift m the strangeness S = +1 channel for rriK = tti^ (left panel) and mx rn-Tr (right 
panel) as functions of the kaon momentum fc^ = o;^ — mj^ in the BSA. 

B. Fluctuations orthogonal to rigid rotations 



Similar to the above study we eventually have to consider small amplitude fluctuations about 
the rotating soliton in the RVA, i.e. rja in eq. (|2.1H) in sections V-VIII. In these sections we will 
show that the phase shifts extracted from rj and ry are identical in the large- A^c" limit; the ultimate 
test on the RVA. In contrast to the BSA fluctuations r/^ explored so far, the rj^ must be constrained 
to the subspace orthogonal to the rigid rotations of the hedgehog to avoid double-counting. These 
constraints follow directly from the ansatz, eq. (|2.11j) and we will derive them in the appendix. 
Here we will study the solutions to the equations of motion that arise when the BSA equations of 
motion are supplemented by these constraints. We call these solutions rj and 11. 

As noted above, the RVA Lagrangian contains terms that are linear in rj providing Yukawa 
couplings to rotational excitations of the soliton. The equation of motion for rj will then differ from 
the one for r/ in two aspects: (i) the constraints and (ii) exchanges of resonances induced by the 
Yukawa couplings. The important point here is obvious: supplementing the equation of motion for 
T] by the constraints or omitting the Yukawa exchange terms from the one for rj identically gives 
the one for the auxiliary field rj. Since r] is obtained by removing the resonance pieces from rj it 
is the background contribution to the scattering amplitude. This strongly motivates the study of 
the auxiliary field fj as an intermediate step to construct the RVA fluctuations, rj; even though the 
introduction of rj seems somewhat artificial at this moment. 

The constraints only concern the P-wave channel and read 

^= r'^dr z{r) \Tl{r) - X{r)r]{r)] = and X= r'^dr z{r)MK{r)r]{r) = . (3.12) 
Jo Jo 

We add these constraints to the Hamiltonian with global Lagrange-multipliers a and (3. From this 
we obtain the homogeneous equations of motion for the fluctuations 11 and rj, that are orthogonal 
to the rotations of the soliton, as modifications to eqs. (|3.3|) 



n(r) = [ujMKir] - X{r)]ri{r) + aMK{r)z{r) 



rj{r) + 



A^(r) 
MK{r) 



A(r) 
MK{r) 



n(r) - aX{r)z{r) - j3MK{r)z{r) . 



(3.13) 



We multiply the above equations by z{r) from the left and integrate over the radial coordinate to 
determine the Lagrange-multipliers using the constraints, eq. ()3.12|) 



a 



CO fOO 

2 J /„\ Inx /^M ^\ „„J a , . „, I t / 2, 



r dr z{r) [2X{r)]r]{r) and /3 = + (m^^ — m^) / r dr z{r)r]{r) , (3.14) 
Jo 
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BOUND STATE APPROACH 
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FIG. 2: The full lines are the same phase shifts as in figure ^ The dashed lines are the background phase 
shifts in the strangeness 5 = +1 channel as functions of the kaon momentum /c^ = cj^ — as computed 
from the integro-differcntial equation H3.15|l . The dotted lines give their difference to be identified as the 
resonance phase shift. Left panel uik = tIt, right panel uik TO,r- 



since J r'^drzMxz = 1 and J r'^drz2Xz = ujq. We eliminate 11 from the second equation ()3.13() . 
substitute the above Lagrange-multipliers and obtain the linear second order integro-differential 
equation for rj, 



h^fi{r) +uj[2X{r) - ujMK{r)]fi{r) 



[2A(r) -{uj + wo) Mx(r)] z{r) / r'^dr' z{r')2X{r')r]{r') 

Jo 



niZ) I r'^dr' z{r')rj{r') . 
'o 



(3.15) 



Any solution to this equation automatically satisfies J r'^drzMxr] = 0. This is easily confirmed by 
multiplying eq. (|3.15|) from the left with z(r) and integrating over the radial coordinate. Obviously 
this integro-differential equation projects eq. ()3.5() on the subspace orthogonal to the rigid rotations 
of the hedgehog. As a result the integro-differential equation does not have a bound state solution 
to be associated with the A hyperon. In figure El the dashed lines are the phase shifts 6 for 
strangeness S = +1 as computed from the integro-differential equation (|3.15|) . This phase shift is 
quite shallow with an indication of a small repulsive potential. As explained at the beginning of 
this subsection, rj parameterizes the background part. That is, S is the background phase shift in 
the strangeness S" = +1 channel. On the other hand the total phase shift, 6 {cf. figure ^) is the 
sum of the resonance phase shift, 6r and the background phase shift. 



SRik) = 6{k) - 5{k) 



(3.16) 



We also display this resonance phase shift in figure 121 Obviously it exhibits a clear resonance 
structure as 5R{k) smoothly passes through 7r/2. Thus we reach a central result of our studies: 
Once the background phase shift is properly identified and subtracted, chiral soliton models indeed 
predict a distinct resonance structure in the strangeness 5 = +1 channel even and in particular 
for Nc oo. In section V we will show that in the limit Nc — > oo it is exactly (5r that arises 
from exchange of the O"*" pentaquark (and the A hyperon). Figure [21 suggests that the resonance 
is quite broad. In section VI we will discuss that finite Nc effects are substantial in sharpening 
this resonance. Note also from figure (21 that the position at which 5R{k) passes tt/2 is somewhat 
lower than predicted by the estimate, eq. ()3.11|) . In section VI we will explain that this pole shift 
originates from the well-know effect of iterated resonance exchanges. Actually this is standard in 
scattering theory and does not alter the observation that loq is the "bare" resonance or pole mass. 
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C. Hyperfine splitting 



Having established the existence of states in both the 5 = — 1 (as bound state) and 5 = +1 (as 
resonance) channels we would like to complete this section discussing the hyperfine splitting in the 
BSA for later comparison with the RRA. Even more, we have seen that the corresponding kaon 
wave-functions are dominated by the zero mode component, z(r). A posteriori this justifies the 
estimates for the mass differences to the nucleon, eqs. (|3.10j) and (|3.11|) . These estimates concern 
the large- A'^c' limit and are thus degenerate in spin and isospin. To lift this degeneracy, states with 
good spin and isospin need to be constructed in the framework of the BSA. For the S* = — 1 channel 
this has been performed in ref. [i^ and we will generalize it to also apply for 5 = +1. Starting 
point is the parameterization, eq. H2.1(J() with the zero mode wave-function substituted for r]a in 
the P~wave, cf. eq. (|3.2)) . Spin and isospin states are then generated by canonical quantization of 
the SU{2) collective coordinates, A2. Substituting this parameterization into the action adds 



dK{t) 
dt 



TK{t) - K\t)f 



dK{t) 
dt 



Nc 



n-K^{t)TK{t) (3.17) 



to the Lagrange function |49|. Here \/2K{t) = J da;x(w)e*'^* contains the amplitudes for the kaonic 
states in the S = ±1 channels. The dependence on the collective coordinates is parameterized by 
the angular velocities 



4(i)^ = if.ji. 

' dt 2 



In addition the moment of inertia for spatial rotations appears. 



= - / dVsin^F 



fl + 



1 



sin^F 



5.115GeV 



Treating this additional term perturbatively modifies the Hamiltonian by 



1 



29^ 



2Qk 



nV A' - ArVn ) + 



AQk 



K^tK 



(3.18) 



(3.19) 



(3.20) 



where Jj = and 11 is conjugate to . The computation of matrix elements of AH is 

standard [i^ and yields the hyperfine contributions to the masses of baryons with strangeness 
S = ±1, 



AMs 



1 



26, 



csJ ( J + 1) + (1 - cs) /(/+!) + -cs {cs - 1) 



(3.21) 



with J and / being the spin and isospin quantum numbers of the considered baryon, respectively. 
The hyperfine parameters 



1 



and 



1 



(3.22) 



80 KUJ A + Nc ^ SOkujb-Nc' 

are 0{N^) so that the coefficient in eq. ()3.21|) ensures limiy^^oo AMs = 0. We remark that 

the expressions in eqs. 1)3. 10() . ()3.11|) and (|3.22|1 are the BSA predictions for fluctuations, r] in the 
subspace of rigid rotations. 

To this end, the mass differences of the A-hyperon and the G"*" pentaquark with respect to the 
nucleon read 

3 



Ma, 



e 



M 



N 



WA,e + AMs 



89, 



(3.23) 
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Similar relations can be easily written down for other S = ±1 baryon states. The last term in 
eq. (|3.2H) is often omitted as it contains terms that are quartic in the fluctuating field and the 
calculation is not completely carried out at this order. The comparison with the RRA suggests 
that these quartic terms contribute 9/8@k to the mass of the 5=1 baryons. This is not negligible 
for Nc = 3 and clearly shows that subleading contributions must be considered for predictions to 
be reliable. These subleading contributions can be evaluated within the RRA that we will discuss 
in the next section. 



IV. RIGID ROTATOR APPROACH (RRA) 

In this section we will explain how baryon states emerge as soliton excitations that are members 
of flavor SU{3) representations. They are obtained from canonical quantization of the collective 
coordinates ^3 in the ansatz, eq. (|2.11() . For this purpose we do not need to consider the fluctua- 
tions, rja and relegate the discussion of the coupling between ^1^ and t/q to the following section. 
The resulting Lagrangian for the collective coordinates reads [30| 

L = -M+ ^e^l + Ibk^I - ^ - ir (1 - Dss) . (4.1) 

The soliton mass, M is the functional which upon minimization determines the chiral angle, F{r). 
Its particular form is of no relevance here because we will only consider mass differences. The 
moments of inertia, Q^r, @k and the symmetry breaking parameter, T have already been introduced 
in the previous section, cf. eqs. (|8.6|) . (|.S.9|) and (|3.19|1 . In equation (|4.1j) also angular velocities 
for rotations into strangeness direction appear. These are straightforward generalizations of the 
definition irTTsll to ^3 

4W^ = ^AA. (4.2) 

Finally the symmetry breaking part of the action, eq. (|2.6j) causes the explicit dependence on the 
collective coordinates. This is conveniently displayed with the help of their adjoint representation, 



1 

-tr 
2 



(4.3) 



We note that in more complicated soliton models (e.g . with vector mesons or chiral quarks) the 
symmetry breaking part of L has a richer structure [SJI which eventually leads to better agreement 
between predicted and observed baryon spectra. However, we presently compare two different 
treatments (BSA vs. RRA) of the same model and thus other symmetry breaking terms are of no 
relevance to us. For quantization we compute the conjugate momenta. 



dL 



— Qx^a , 
Nc 
2V3 ' 



-Jr 



a 
a 
a 



« = 1,2,3 



a 



4 7 



(4.4) 



The quantization prescription then demands the commutation relations [Ra,Rb] = —ifabcRc with 
fabc being the antisymmetric structure constants of SU{3). That is, the Ra are the generators of 
the intrinsic SU{3) group. The hedgehog structure of the classical soliton ensures that Jj = —Ri, 
{i = 1, 2, 3) is the total spin operator. The Wess-Zumino term, eq. (|2.5|) plays a decisive role. It is 
non-local and only linear in the time-derivative such that it contributes the term linear in to the 
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collective coordinate Lagrangian. Therefore the corresponding conjugate momentum is constrained 
to Rg = and we require a Lagrange multiplier, fj, for this constraint when constructing the 
Hamiltonian for the collective coordinates by Legendre transformation from eq. ()4.1|) 



H = M + + -L-Rl + ir(l - - , (^R, - ^) . (4.5) 

Demanding that L can be obtained from H by the inverse Legendre transformation requires to 
set fj, = Qs- Since [H,Rs] = this angular velocity remains undetermined and different choices 
merely imply different gauges. This vanishing commutator also suggests that the constraint can 
be implemented on the eigenstates of H. Some of the discussion on finding these eigenstates for 
arbitrary Nc is already contained in the literature ji^. Here we focus on the comparison between 
BSA and RRA for the S = +1 pentaquarks, which is new. So is the non-perturbative investigation 
of symmetry breaking effects in the RRA for arbitrary Nc^. 

We concentrate on the case of vanishing symmetry breaking, F = for a first discussion on 
finding the eigenstates of the collective coordinate Hamiltonian, eq. 1)4. 5|) . For F = the eigenstates 
of H are members of definite irreducible representations of SU{3) that are characterized by the 
labels {p,q)- The eigenvalues of H can then be computed from the eigenvalue of the quadratic 
Casimir operator of SU{3), 

The constraint on Rg only allows representations with 

p + 2q = Nc + 3t (4.7) 

where t = 0, 1, 2, . . . refers to Biedenharn's triality 0]. The physical hypercharge for baryons is 
generalized to y = ^ + S. Note that because of l^ax = ^ + t this triality equals the largest 
strangeness quantum number ^max = t contained in the representation {p,q)- The eigenstates of 
H are most conveniently arranged according to increasing t. For t = they are 

ip,q)=(2J,^^^) (4.8) 



where J = |, |, . . . , ^ are the allowed spin eigenvalues, i.e. (Rf) = J{J + 1). From equation (|4.()|1 
we compute the corresponding energies 



which coincides with the SU{2) mass formula [lj| up to an additive constant for nucleon and A 
type states. For i = 1 we find 



^ A simplified treatment of symmetry breaking to only leading order in perturbation theory gives results consistent 
with ours Is^ . 
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with J = i, |, . . . , ^ + 1. In particular the lowest lying multiplets with spin J = h and J = ^ are 



J=- : ip,q)=i3, 2 



2 

iVc-3 



2 

iVc + 1 



'8", "10",... 
10", "27",.. 



(4.11) 



where the expressions to the furthest right refer to the dimensionalities of the representations for 
Nc = 3. The above assignments are unique in soliton models and so are the flavor quantum 
numbers attributed to the lowest lying baryon states within the multiplets. Essentially this is 
inferred from the fact that these assignments minimize the energy functional. For example, for 
large- A'^C' the symmetry breaking term contributes 



'1-D. 



88, 



r 

2 



1 



Nc' 



-Y 



3T 

'2N~c 



S 



iS<t) 



(4.12) 



to the energy functional if symmetry breaking was treated in first order perturbation. For finite 
Nc this generalizes to the fact that symmetry breaking causes the masses of baryons within a 
multiplet to decrease with increasing strangeness. For t = we therefore obtain the lightest J = i 
baryon, the "nucleon" with zero strangeness and isospin 1=2 and the lightest J = 2 baryon, the 



" A" with zero strangeness and isospin I 
and J 



For t = 1 the lowest lying "pentaquarks" with J 



2 carry strangeness S* = +1 and isospin 7 = and 1=1, respectively. Of course, for 
finite symmetry breaking mixing with higher dimensional multiplets occurs (see below). Although 
such non-linear effects may not be omitted, the spin-flavor assignments are unaffected. 

The lowest lying pentaquarks with S = +1 dwell in t = 1 representations and the two mass 

HI 



formulas (|4.1flj) refer to the spin-isospin relations / 
corresponding energies 



M + 



J(J + 1) 
29^ 



+ 



+ 



1 



Nc 

2Qk ' ^Qk 



1(1 + I) 



It is straightforward to derive the 



(4.13) 



The terms 0{1/Nc) exactly match^ those in eq. (|3.21j) when equating the hyperfine energies that 



are related to expressions quartic in the fluctuations to 



Even more, we may compute the 



mass difference between the nucleon (in "8") and pentaquarks with 7 = and J = o "10 



E-. 



10 



Eh 



iVc + 3 



(4.14) 



which coincides with ujq as Nc -^00, as it should. Thus we conclude that the BSA and RRA are 
consistent when flavor symmetry breaking is omitted. Note that this mass difference acquires a 
factor 2 for iVc = 3. 

We will next consider the physical case of non-zero flavor symmetry breaking. As a flrst measure 
a perturbative treatment of flavor symmetry breaking seems temptin g be cause it reproduces the 
Gell-Mann-Okubo mass relations p^ . IH^ in the RRA for Nc = 3 |3ri|| . However, we already 
argued in the context of eqs. 1)3. 1U() and ()3.11|) that the perturbation series in F does not converge 
for the BSA. We recall that many of the pentaquark studies in the RRA treated flavor symmetry 
breaking perturbatively (at best to second order) 0, 0]®. This seems at odds with the BSA. 



Without symmetry breaking we have lob ~ Nc/'^Qk and thus 1 — c+ = QtJQk- 

In the case of pentaquarks symmetry breaking was treated exactly in refs. |flL 1^ l38l l59| . 
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mass differences, 0{fi^), S = - \ 



mass differences, O(N^), S = +l 
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e - N (RRA) 
a - N (BSA) 



FIG. 3: Mass differences at 0{N^) computed for e = 4.25 in the Skyrme model as functions of Nc- In 
the RRA they are the corresponding differences of the eigenvalues in eq. H4.15|l while in the BSA they are 
extracted from eq. H3.23|l . Left panel AS = — 1; right panel AS' = +1. 



Also, the Gell-Mann-Okubo relations are not reproduced in the BSA ISII 
the eigenvalue problem 



Hence we have to solve 



M + 



J(J + 1) 
26^ 



+ 



2 



1 - Dss) >^ = together with Rs^ 



Nc 
2^3 



(4.15) 



(numerically) exactly. For Nc = 3 this has been accomplished some time ag o l6Ci a nd frequently 
used since then for the computation of a variety of baryon properties, cf. refs. |59l.l6ll. l62ll . This can 
straightforwardly be gener alized by putting Yr = 2Rs/V3 = Nc/3 in eq. (2.18) of ref. |60|, together 
with the Nc dependencies of G^, cf. eq. 1)3. 6() and T (paragraph after eq. p.lUI) '). Equivalently, 
the operator Dgs may be diagonalized in the space of SU{3) representations with i?8^ = 
using SU{3) Clebsch-Gordan coefficients. The only condition on Nc is that it must be taken odd. 
In section VII and VIII we will employ the resulting wave-functions to compute matrix elements 
of collective coordinate operators for arbitrary odd Nc and non-zero symmetry breaking. 

For the previously mentioned parameters we present the resulting mass differences in figures El 
and0I In figure|31we concentrate on mass differences that scale like 0{N^), i.e. between baryons 
whose strangeness quantum numbers differ by one unit. As Nc — > oo the BSA for fluctuations 
in the rotational subspace predicts uja and ujq for those mass differences regardless of spin and 
isospin. Obviously this is matched by the RRA. In figure |1] we display the mass differences that 
scale like 0{1/Nc), i.e. we compute the hyperfine splitting in both approaches. In doing so we only 
consider combinations for which in the BSA the ordering ambiguities and omissions from 0{rj'^) 
terms cancels, i.e. baryons with identical strangeness. Again, perfect agreement between the two 
approaches is observed for Nc — > oo. However, at Nc = 3 sizable 1/Nc corrections occur that 
cannot be accounted for by the BSA. 

A couple of further comments are in order. First, when computing matrix elements we observe 
that, cf. eq. (jHE]), 



0{N[ 



(4.16) 



that is, the net effect of flavor symmetry breaking in the RRA is an 0{N^) contribut ion to the 
baryon mass. Again, this is consistent with the BSA where flavor symmetry breaking enters the 
0{Nq) equation of motion for the kaon fluctuations (|3.5|1 . Second, the literature contains claims 
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mass differences, 0(1/N„) 



mass differences, C>(l/N„), S = -l 



IBSg88a8SSS«( 



A - N (RRA) 
A - N (BSA) 
0*- (RRA) 
0*- (BSA) 




FIG. 4: Mass differences at 0{\/Nc) computed for e = 4.25 in the Skyrmc model as functions of Nc- Left 
panel: baryons with strangeness = 0, 1; right panel 5 = — 1. Sec also caption of figure |3| 



that the A-nucleon mass difference would scale like 0{N^) once symmetry breaking is incorporated 
in the RRA |lll |. The actual calculation, displayed in figl3 shows that this is not correct. As in 
two extreme cases, zero and infinite^ symmetry breaking, this mass difference scales like C'(l/A''c). 
The physical value, rriK = 495MeV is intermediate to these extreme cases and it thus comes to no 
surprise that this mass difference vanishes as Nc — > oo. 

The second remark concerns the eventual restriction of collective coordinates to only dynamical 
zero modes. For zero symmetry breaking (F = 0) the rigid rotation eigenstates of ea. H4.15|) are 
states in the 8, 10, 10, 27,... dimensional representations of flavor S'C/(3). Though the anti-kaon 
fluctuations are zero frequency solutions in the BSA, the kaon fluctuations are not due to the 
Wess-Zumino term, i.e. the latter are not really dynamical zero-modes. Therefore it has been 
suggested to introduce collective coordinates only for the anti-kaonic zero-modes in a manner that 
only the octet and decuplet would be treated collectively [31 • Such an approach seems complicated 
because constraints on both the fluctuations and the collective coordinates are required. Putting 
the question of feasibility aside, there is no advantage in doing so to study the pentaquark since 
such a treatment essentially corresponds to the BSA and prevents the incorporation of finite Nq 
effects. The restriction to the octet and decuplet seems even more problematic in view of flavor 
symmetry breaking (F 7^ 0) as it mixes baryon states from different SU{^) representations |6ll ]. 
In particular the nucleon wave-function contains admixture of the nucleon type state in the 10. 
Thus the latter may not be omitted. 

To round up this section on the RRA we note that other approaches to quantize chiral solitons 
in flavor 5C/(3), like the slow rotator [^^l or the breathing mode [H^l approaches, are essentially 
extensions of the RRA. 

So far we have established the relation between rotational excitations of the soliton and the 
bound states in the \aige-Nc limit. In order to study the decay width of the pentaquark we 
need to introduce fluctuations by extending the RRA to the RVA. 



® This is the SU (2) version of the chiral soliton model. 
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V. COLLECTIVE ROTATIONS AND VIBRATIONS 




We will now include the kaon fluctuations in the RVA according to the ansatz (|2.1H) and 

compare it to results for the fluctuations in the BSA, section III. The flavor rotating hedgehog 
A3{t)Uo{x)A^^{t) is not a solution to the field equations. Hence the expansion of the action with 
respect to rj has a linear term. This linear term gives rise to Yukawa couplings like KNA and 
KNQ. As a consequence the KN phase shift in the isospin 1 = channel^'' has contributions from 
(iterated) exchange of A and Q. Technically, these exchange contributions emerge in the equation 
of motion for 77 as a separable potential of the form " (Yukawa interaction) x (propagator) x (Yukawa 
interaction)". The reader may want to peak ahead to eq. for the explicit expression. We 

have already noted that the background field rj, that we have computed in section III, obeys an 
equation of motion that is obtained from the one for rj by omitting the Yukawa couplings. That is, 
the respective equations of motion differ by this separable potential and we can straightforwardly 
construct rj from rj using standard Lippmann-Schwinger equation techniques. This will make the 
role of exchanged resonances very transparent. For convenience we list the notation for the involved 
radial functions 

r]{r) : solution to eq. ()3.5|1 in the BSA 

r]{r) : solution to the integro-differential equation, 1)3. 15(1 with constraints 

rj{r) : solution to the integro-differential equation, l|5.1Uj) with constraints and Yukawa exchange 
We will establish two central results in the large-A'^c limit: 

A) We compute the RVA fluctuations rj by supplementing the equations of motion for r] with 
the effects due to A and G exchanges. We compute the phase shifts for rj and show that they 
are identical to those computed from rj in the BSA, eq. (|2.1()|) . 

B) We separate the Yukawa exchange contributions to the scattering amplitude and find the 
corresponding phase shift to equal the resonance phase shift displayed in figure [21 

In doing so we need to identify the Yukawa couplings by computing the terms linear in rj. Note that 
in the action all linear terms must occur as iso-singlets. Then there are two sources for Yukawa 
terms: 1) the flavor symmetric piece Tsk + ^wz has couplings to the angular velocities $7^, 2) the 
symmetry breaking piece has a contraction between Dga and rja- For the moment we will only 
consider the first case and relegate the discussion of flavor symmetry breaking to section VII. 

In the flavor symmetric case there are already many rotation-vibration coupling terms linear 
and quadratic in $7^. In order to keep the calculation feasible we keep only those terms that 
survive in the large- A'^c limit. The so-selected Yukawa coupling is then treated to all orders 
in Nc- Admittedly, this procedure is not completely consistent and the neglected couplings are 
not necessarily small at Nc = 3. However, a major purpose of the present investigation is to 
show the equivalence with the BSA. For this purpose it completely suffices to only consider the 
leading Nc couplings. In order to isolate these terms we need to discuss the Nc scaling of the 
angular velocities 0,a to identify those Yukawa terms that survive in the large- A'^c' limit. That is 
straightforward for a = 1, . . . , 7 as we simply invert eq. (|4.4|1 . 




ni = 0{l/Nc) I.e. R^ = 0{N^) 



The I — 1 channel, which we discuss in the last paragraph of section VIII, has E* and O* exchange contributions. 
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E 



C2 - J(J+1) - 



= ^ ^^^2^ ^ ^ 0, = O ( 1/ ViVc ) i-e. Ra = 0{VNc), (5.1) 

K 



as the moments of inertia ©tt and 0^- are O (Nc)- In the large- A'^c' expansion (i = 1,2,3) is 
thus subleading to (a = 4, . . . , 7) and may be omitted. As shown in the appendix, the angular 
velocity may be absorbed in a suitable redefinition of the energy scale 

a/3 

co^u + ^ns, (5.2) 

where lv is the frequency of the kaon fluctuation, cf. eq. ()3.2|) . Thus, to compute the O-r/ Yukawa 
coupling we only need to consider the kaonic angular velocities = —Ra/&K with a = 4, . . . , 7. 
This lengthy calculation described in the appendix not only requires to extract the terms in the 
action that are linear in r/ but also to consider the constraints, eq. (|3.12|) because they contain pieces 
that are of the same linear order in rj. Having extracted these linear terms we are unfortunately 
not yet in the position to compute matrix elements of the interaction Hamiltonian (|5.5|) between 
and KN states in the space of the collective coordinates. The reason is that rj parameterizes 
the kaon field in the intrinsic frame. A collective flavor rotation relates the intrinsic fluctuations 
to those in the laboratory frame, that are observed in kaon-nucleon scattering, 

= Dabrjb , (5.3) 

with o, 6 = 1, . . . , 8. In the large-A'^c' limit, D^i and Das, the couplings to the intrinsic pions and 
r/, respectively, are suppressed. Hence we approximate 

Ca = Daf^rjfS , (5.4) 

with a, /3 = 4, . . . , 7. Finally the resulting interaction Hamiltonian reads. 

Hint = -pr^ dia/B D^aRfB d^r sm—— [2X{r) - uJoMK{r)] Xi^^{x ,t)Ri3 , (5.5) 
fx J 2 

where diap are the symmetric structure constants of SU{3). The pieces involving loq = Nq/^Qk 
and A(r) stem from Tsk and Twz, respectively. From eqs. (|2.7|) . (|5.H) we easily verify that 

^int = O(iV^). 

We need to compute matrix elements of the interaction Hamiltonian, eq. (|5.5j) between a kaon 
nucleon state coupled to good spin (J) and isospin (/) and a collective excitation with the same 
quantum numbers. We denote the operator that annihilates a kaon in the laboratory frame with 
isospin projection i = zt^ and orbital angular momentum {I, m) by a^*^ and use the Wigner-Eckart 
theorem to reduce the matrix elements of i^int to those computed in the space of the collective 
coordinates, 

- 2 

where 7(±|) maps the isospin projection onto the flavor indices, D±a = i iD^a- The arrow 
indicates the spin projection J3 = i. In what follows we will always adopt this value without 
explicit mention. Orbital angular momentum £ = 1 is projected onto by Xi in Hi^t- Of course, 
that is nothing but the statement that only the P-wave channel is affected by the rigid rotations. 
A similar reduction formula holds for the matrix element of i^int between final A and initial anti- 
kaon nucleon states. In the flavor symmetric case the relevant matrix elements of the collective 
coordinate operators are 



^ {@K!^dmal3D^ii)aRf^\{KN)j^^^j^i) = ^6 (6+ T \dsafBD+^Rp\n T) (5.6) 



in+u n R I \ + 1 ^c + l ^ Nc 

(6 \d,^pD^^Rp\n) = -^^2iNc + 3){Nc + 7) 1^ 
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{k\d^^pD.^Rp\p) 



iVc + 5 



32 



2 {Nc + 7) ^2Nc + 6 

For the matrix elements in eq. (|5.7|) we have factored the large-A'^c value for later convenience 



(5.7) 



/ {Nc + 1)^ 
/ Nc {Nc + 3) {Nc + 7) 

2 {Nc + 5) ^ 
'{Nc + 7) ^Nc{Nc + ?>) 



I 3^5 



Nc oo 
Nc = 3 

Nc^ 
Nc = 3 



oo 



(5.8) 



Only the Yukawa coupling survives as Nc oo. Hence there is only one intermediate state in 
that limit. For the physical value, Nc = 3, the collective coordinate part of the A Yukawa coupling 
constant is actually larger in magnitude than that of 0^. 

Obviously this explicit derivation of the interaction Hamiltonian yields only a single SU{3) 
structure (i.e. di nnXiD nnRpt) for the Q — N transition. This is to be contrasted with earlier 
studies Ha, H 



42| | who "invented" three structures, even in the flavor symmetric case. 
We want to add the particle exchanges induced by H\yit to the differential equation, (|3.15|) for the 
constrained fluctuations, r/ (still omitting symmetry breaking). The propagators are {loq—u})~^ and 
{oj\ + a;)^^ for O and A, respectively, because the differential equations for kaons and anti-kaons 
are related by w ^ — w. Thus the Yukawa terms induce the separable potential 



Ke|//int|(A'iV),=o)|' , mH,,,,\{KN) 
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(5.9) 



Note that these matrix elements concern the T-matrix elements in the laboratory frame. However, 
we observe from eq. (|3.ip that for the channel the laboratory and intrinsic T-matrix elements 
are identical. Hence we may add the exchange potential, eq. 1)5. 9() in the intrinsic frame. We then 
end up with the modified integro-differential equation 



h'^r]{r) + uj [2A(r) - u;Mx(r)] ??(r) = -z{r) 



2A(r) - {u; + uJo)MK{r)-iOo 



r' dr' z{r')2X{r')r]{r') 



Nl 



UJQ — UJ 



UJ 



{2\{r) - uJoMK{r)) 



(5.10) 



in the flavor symmetric case. The particular coefficient, ujq of the Yukawa contribution will become 
clear when we discuss the width of O"*" in the next section. The Yukawa contribution to the 
integro-differential equation is orthogonal to the zero-mode. 



r^dr z{r) [2A(r) — u;oAf/^(r)] z{r) = 



(5.11) 



and thus any solution to eq. (|5.1()j) still satisfies the second constraint in eq. (|3.12|) which we 
subsequently employ to simplify the integral. 

The essential observation is that eq. (|5.10|) has a simple solution as Nc — > oo (Xe = 1 and 
^A = 0) 



rj{r) = r]{r) — az{r) 



with 



drr'^ z{r)MK{r)r]{r) . 



(5.12) 



Here r]{r) is the solution to the unconstrained equation (|3.5j) in the BSA. The radial function, z{r) 
associated with the zero-mode is localized in space. Thus the phase shifts of rj and rj, which are 
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extracted from the respective asymptotic behaviors are identical. This proves om' assertion A). In 
turn this shows that the difference of the phase shifts defined in eq. (|.S.16|) is directly related to the 
resonance exchange. This proves our assertion B). In section VII we will show that these assertions 
also hold when symmetry breaking is included even though that case is significantly more involved 
as e.g. hmTv^^oo -'^A / 0. 

VI. THE WIDTH OF THE 9+ 

In this section we will show that the Yukawa terms can be imposed to compute the width of 
the exchanged particles. For simplicity we will still omit flavor symmetry breaking. 



A. Large A^c limit 



In a flrst step we show that in the large- A^c* limit the iterated exchange, that is induced by 
the Yukawa terms in eq. (|5.10j) . leads to the resonance phase shift discussed already in section IV. 
That is, we want to prove assertion B) by showing that the resonance width associated with 6R{k) 
equals that induced by the Yukawa coupling of 0"*". 

For large-Ac we may omit the A~exchange. Then the reaction matrix formalism is most 
suitable to compute Sii{k) from the iterated exchange of a single resonance. The (real) reaction 
matrix obeys the Lippmann-Schwinger equation 



R = V + VGqR , 



(6.1) 



where Go and V are the Green's function for the constrained problem, eq. (|.S.15|) and the separable 
potential due to the particle exchange, eq. (|5.9() , respectively. We stress that this is an exact solution 
to the scattering problem rather than any type of Born or distorted wave Born approximation. The 
solution to the Lippmann-Schwinger equation is most conveniently obtained in a basis independent 
formal language 



iruj 



q^dq 
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UJ 
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with 



\V-coJiV- 
OJ + OJa 



(r|5e) = Ae (2A(r) - cuoM^(r)) z{r) 



(6.2) 



where ujn 



+ m\ and V denotes the principle value prescription. In the above expression 77 
appears, rather than r/, because the interaction Hamiltonian is treated as an additional interaction 
for the constrained fluctuations and we thus require the Green's function for rj. The potential is 
obtained from eq. (|5.9|) by setting Aa = 0. We flnd the reaction matrix 

^ ^ t^oAgl (2A - ujqMk) z){z (2A - uj^Mk) \ 
u;-u;e-^o(^(2A)|Go(a;)|(2A)z) 

since {r]J\MKz) = 0. The resonance phase shift induced by V is obtained from the diagonal matrix 
element of the R-matrix, k{r]^^J\R\ri^^^) = — tan ((5r(A:)), yielding 



tan (5r(A:)) 



This phase shift exhibits the canonical resonance structure with the width 

2 







r dr z{r)2X{r)r]^^{r] 



(6.4) 



(6.5) 
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and the pole shift 

AK) = - V / qdq 

Of course, we have numerically verified that in the large- A'^c limit with Xq = 1, the phase shift 
from eq. (|().4)1 is identical to the one that we calculated as the difference in eq. (|I-{.1()|) . For finite Nc 
the width F acquires a factor 7^ 1 and the R-matrix formalism becomes two-dimensional (A 
and exchange). 

We stress that although the singularity of the scattering amplitude is shifted by the amount A 
due to iterated resonance exchange, the "bare" mass of the resonance is ujq, the pentaquark mass 
predicted in the RRA. 

Alternatively we have applied Fermi's golden rule to compute the transition rate Q — > KN from 
the interaction Hamiltonian, eq. ()5.5() . This leads identically to the result displayed in eq. ()6.5|) 
and shows that the way we have included the separable potential in eq. (|5.1flj) is unique. 



+ 



■UJq CJfc + UJq 



(6.6) 



B. Finite Nc 

As already mentioned, the R-matrix formalism is two dimensional at finite Nc because the A 
exchange cannot be omitted. It is therefore more convenient to compute the resonance phase shift 
in a fashion similar to eq. (|3.16|) . First we compute the phase shift with both Yukawa interactions 
included from rj in eq. (|5.10j) . From that we get the resonance phase shift by subtracting the 
background phase shift that we extract from in eq. 1)3. 15() . The resulting resonance phase shifts 
in the channel are displayed in figure [3 It is illuminating to see how the BSA result emerges 
in the limit Nc 00. Obviously this resonance becomes sharper and the pole shift decreases (in 
magnitude) to A = — 14MeV as the number of color assumes its physical value Nc = 3. This 
is mainly due to the reduction of Xq. Furthermore, at finite Nc the bare resonance position is 
increased to ujq = {Nc + 3)/4G/^, according to eq. ()4.14p . This A'^c'~dependence of the resonance 
position is clearly exhibited in figure El 

There are two competing effects on the width when approaching the physical value, Nc = 3. The 
position of the pole is increased by about a factor of two. This increases the available phase space 
and thus the width. On the other hand, the collective coordinate matrix element Xq decreases by 
approximately a factor 3, cf. eq. (|5.8j) . According to eq. (|6.5j) this decreases the width. 

At this point it is also illuminating to consider the large- A'^c' expansion of these matrix elements 
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in more detail. For example 

. . 7 52 372 2628 ,^ 

does not converge for Nc < 7. This clearly demonstrates that taking only leading pieces from a 
large- A'^C' expansion of these matrix elements yields incorrect results for Nq = 3. It is unavoidable 
to include finite Nc effects. Any approach |2J] that employs 1/Nc expansion methods for exotic 
baryon matrix elements at Nc = 3 seems questionable. 



C. Comparison with other approaches to compute the width 



We stress that the appearance of only a single rather than two SU{3) structures in the transition 
Hamiltonian in the flavor symmetric limit is not an artifact in the (eventually oversimplified) 
Skyrme model whose action is a functional of only pseudoscalar mesons. Rather it generalizes to 
all chiral soliton models. The reason is simple. Starting point in all these models is the introduction 
of intrinsic fluctuations about the classical soliton, i.e. the generalization of the ansatz, eq. ()2.11|1 
to parameterize additional fields. In this intrinsic formulation the interaction Hamiltonian must 
be an isoscalar operator that is linear in kaonic fluctuations and projects onto the P-wave channel. 
At 0{N^) only diapXirja^p is possible. Of course, in more extended models rja also involves other 
strange degrees of freedom, e.g. K* . In the flavor symmetric formulation the adjoint representation 
of the collective coordinates, Dab, emerges only by the transformation to the laboratory system 
as in eq. (|5.3j) . To see whether or not the limitation to only a single SU{3) structure in the 
transition Hamiltonian is an artifact of the Skyrme model it is instructive to recall how the authors 
of refs. [ifll. I41I1 came to the improper use of several SU{3) structures regardless of flavor symmetry 
breaking. Starting point of those studies is the spatial part of the axial current and its identification 
with the kaon field operator via PCAC. This axial current is commonly computed from the rigidly 
rotating classical field only and in general allows for three SU (3) structures in the fiavor symmetric 
case. One of those structures is directly related to the presumably small singlet current matrix 
element and is usually ignored in the context of the O"*" width. In the simple Skyrme model it 
vanishes exactly. However, the other two structures are also present in the Skyrme model axial 
current [g^. They simply do not show up in the transition Hamiltonian. Stated otherwise, the axial 
current operator is not suitable to compute hadronic decay widths in chiral soliton models. A simple 
reason is that the isospin embedding of the hedgehog soliton, eq. (|2.8)) breaks fiavor symmetry. In 
particular this prevents using SU{3) to relate various hadron decay widths among each other. 
For example, the computation of the A ttN width requires a term linear in the intrinsic pion 
fiuctuations which involves collective coordinate operators that are not related to diapXifja^p. 

Of course, in our approach d^A^^^ = fT^m'^cf)"' + 0{4>'^), i.e. PCAC, holds as well. It is nothing 
but the equation of motion for the chiral angle and the meson fiuctuations that we solve. That 
is, in our approach PCAC is satisfied at 0{Nc) and 0{N^) while in refs. [ifll. l4ll| the equation of 
motion is only considered at 0{Nc). The main point, however, is that within soliton models the 
right hand side may not be the interpolating field for the pseudoscalar meson in the final state of 
hadronic decays. 



VII. SU(3) SYMMETRY BREAKING 



In section V we have established the equivalence of the bound state and rotation-vibration 
coupling approaches for computing the kaon nucleon phase shift in the large- A'^c' limit in the flavor 
symmetric case. Here we will extend that discussion to the symmetry breaking case. 
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The interaction Hamiltonian that describes the Yukawa exchange for fluctuations in the labo- 
ratory frame, cf. eq. (|5.5j) . acquires an additional term from the symmetry breaking part of the 
action, Fsb; 



int 



@K 



dia/sD-yaR/S 



F(r) 

dVsin— ^ [2A(r) 



4/. 
V3 



diapD-yaDsiS 



^3 ■ Fir) 
(Tr sm^^ 



sr 

89^ 



MK{r) 



C-y{x,t)Xi. (7.1) 



We provide the detailed derivation of this interaction Hamiltonian in the appendix. In principle, 
there are many pole contributions, such as from states with 5 = ±1 and 1 = that stem from 
higher dimensional SU{3) representations. They lead to sharp resonances at higher energies and 
may be omitted for the energy regime of current interest. We thus keep only the contributions 
associated with the exchanges of A and G"*". Then the separable Yukawa exchange potential in the 
P-wave has the general form 



{r\gu) 



UJQ 



\90){g' 



'el 
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sr 
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MK{r)]z{r), u = Q,K. (7.2) 



Note that {z\gy) = 0. This implies that the solutions to the constrained equation of motion (|3.15)1 
with this potential added, also do satisfy the constraint {z\MK'f]) = 0. 

The matrix elements Xjy have already been defined in eq. (|5.7j) . Similarly we introduce 



8Nc 



8Nc 



(7.3) 



where Yq and Y\ tend to unity in the flavor symmetric case as Nq ^ oo. Of course, for non-zero 
symmetry breaking the matrix elements and Yi, have to be computed numerically and cannot 
be given as simple functions of Nc, results are listed in table Q] 

The integro-differential equation with the Yukawa pieces included is a bit more complicated 
than in the flavor symmetric case, eq. (|5.1()|1 . 
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(7.4) 



For simplicity we have omitted the arguments of the radial functions. 

Let rj{r) again be a solution to the BSA equation 1)3. 5() . In the limit Nq oo we require that 
rj = r] — az with a = drr'^ z{r)MK{r)r]{r) is a solution to eq. (|7.4j) also in the case that symmetry 
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breaking {mx 7^ JtItt) is included. Modulo phase conventions this imphes 
Xq = —Yq , X\ 



UJA 



UJQ 



1a. 



Ye = YA 



UJQ +UJA 



(7.5) 



as Nc — > 00. And indeed, the matrix elements in eqs. 1)5. 7() and (|7.3|) give these results when 
computed with the exact solutions to eigenvalue equation (|4.15|) with flavor symmetry breaking 
included! In table Q] we display numerical results for these matrix elements for the physical case 
Nc = 3 and the large- A'^c limit for various values of uik- The third entry, mx = 750MeV yields 
the value of the symmetry breaking parameter F that we would have obtained if we had included 
symmetry breaking terms that contain derivatives of the chiral field. Obviously, there are significant 
deviations from the large- A'^c limit at A'c = 3. 



Nc 


Xe Ye 


Xa Ya 


Xe* Ye* 


rriK = m-K 3 
00 


0.596 0.149 
1 1 


-0.377 0.189 
1 


0.257 0.138 
0.667 0.667 


mK = 495MeV 3 

00 


0.690 0.168 
1.069 0.691 


-0.427 0.195 
-0.378 0.691 


0.307 0.142 
0.713 0.461 


mK = 750MeV 3 

00 


0.801 0.188 
1.157 0.575 


-0.497 0.203 
-0.582 0.575 


0.376 0.147 
0.771 0.383 



TABLE I: Matrix elements of collective coordinate operators for Nc = 3 and Nc 00 and various values 
of rriK, cf. eqs. (|5.7|l . (|7.3|) and (jS.fifl . The third entry, niK — 750MeV is chosen such as to compensate for 
the omitted derivative type symmetry breaker. 

The large-A'c' equation of motion for rj has a bound state solution in the anti-kaon channel 
exactly at the BSA energy uja = 196MeV (for mx = 495MeV). This is a starting point to 
discuss corrections to the RRA spectrum caused by the rotation-vibration coupling for finite Nn 
which improves on earlier attempts that introduced fluctuations induced by the rigid rotations |47|. 
However, that discussion is beyond the scope of this paper. 



VIII. SOLITON MODEL PREDICTIONS FOR 0+ 

We have already noted that in the flavor symmetric case the excitation energy of is altered 
according to eq. ()4.14|) for finite Nc- The excitation energy of A remains at zero. For non-zero 
symmetry breaking the situation is more complicated as both excitation energies depend on the 
value for uik- That is, we need to substitute 

UJA — Ea — En and use — > Ee — Ej\[ , (8.1) 

where Ei\[^e,A are the solutions to the eigenvalue problem, eq. 1)4. 15() in the respective channel for a 
prescribed Nc < 00. In addition, we employ the corresponding eigenstates of eq. (|4.15j) to evaluate 
the matrix elements in eqs. ()5.7|) and 1)7. 3() to numerically compute Xe, XA,Ye and Ya at any value 
of A'c and mx , cf. table We then compute the resonance phase shift as described in section 
VI B. The results are shown in figure El Again these phase shifts exhibit a pronounced resonance 
structure. As in the fiavor symmetric case, finite Nc effects cause the pole position to move to 
higher momenta and also the resonance to become sharper. Although the phase shifts displayed in 
figures El and as a function of momentum seem to be similar this is not the case for the resonance 

energy since it is subject to the dispersion relation u = Jk? + m?^ and rriK differs significantly 



25 




kaon momentum [GeV/c] 



FIG. 6: The resonance phase shift as a function of Nc for rriK = 495MeV. 

in the two case. From figure IHl we read off a G-nucleon mass difference of about 750MeV which 
is a bit higher than the empirical value of about 600MeV if the 0"*" signal around 1540MeV were 
indeed established. We remark, however, that the actual prediction for this mass difference is quite 
model dependent, e.g. we have omitted effects due to fx / f-w- Also, the inclusion of other fields 
like scalar- or vector mesons or chiral quarks will alter the quantitative results. 

The width function, T{u)) is computed from generalizing eq. (|6.5j) to also contain symmetry 
breaking effects, 



' '-Wo 



I) / r^drz{r)ri^^{r) 



i.2) 



which arise from the 0"*" matrix elements of the separable potential, eq. (|7.2j) . Again, r] rather 
than rj appears in the above expression. Note that contributions of the form J drr^z{r)MK{r)fj{r) 
vanish due to the second constraint in eq. (|3.12j) . 

We have already mentioned that chiral soliton models only provide qualitative insight in the 
baryon mass spectrum. At best, we expect the model predictions for masses of exotic baryons 
to be reliable at the, say, lOOMeV level. Let us therefore now assume that the 0^ resonance 
indeed corresponds to the recently asserted pentaquark of mass Mq = 1540MeV and estimate 
its width from r(ti;). First, we have to find the corresponding kaon momentum. This is not 
without ambiguity because in the soliton picture baryon masses are of higher order in Nc than the 
meson fluctuations and thus considered infinitely heavy. Therefore recoil effects are not necessarily 
included. This yields k = 340MeV for the kaon momentum. Of course, a simple calculation using 
relativistic kinematics incorporates recoil effects and leads to k = 270MeV. For these two momenta 
we respectively read off widths of 47MeV and 40MeV for niK = 495MeV from figure[71 Fortunately 
the dependence of r(ti;) on the effective strength of symmetry breaking seems only moderate in the 
resonance region. Hence the extracted width is not too sensitive on the uncertainties stemming 
from the omission of the derivative type symmetry breaker and we finally estimate F ~ 40MeV as 
the width of the 0"*" pentaquark. Though we do not expect chiral soliton model predictions for 
the widths to be reliable at the lOMeV level, we note that this result should be considered a small 
number as it has to be compared to widths of other hadronic decays of baryon resonances, e.g. 
F(A ttN) « 120MeV, empirically. 

Finally let us briefly comment on the width of the 0* pentaquark, the ^ = | and 1 = 1 partner 
of O"*". In large- A'^c the starting point is again the intrinsic T-matrix of the BSA. However, for 0* 
the recoupling between intrinsic laboratory frames is more complicated than eq. (|3.1jl . There are 
two intrinsic T-matrix elements Tee = Ti, and Ts, with different grand spins which contribute to 
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FIG. 7: Model pre diction for the width, r{uj) of 6+ (left) and 9*+ (right) for TVc = 3 as function of the 
momentum k = \J l^P- — nij^ for three values of the kaon mass. Note that the two figures have different 
scales. 



the scattering in the laboratory frame 2^, 5^ . 



i.3) 



While the collective excitation and thus the relevant phase shift displayed in figures ^ and El emerges 
in the intrinsic Q = \ channel, the G = ^ contributes only to the background phase shift and is 
irrelevant for the G* resonance structure. From eq. (|8.3() we notice a factor | not present in the 
corresponding relation for eq. (I.S.lj) . 

In the notation of the RRA and for mx = m.,^ the 0* dwells in the "27"-plet of flavor SU{3) 
for finite Nc while for Nq ^ oo it becomes degenerate with 0+. For Nq = 3 the Q* pentaquark 
is expected to be roug hly lOOMeV heavier than 9+ |ii,Ei- The Yukawa interactions of the RVA, 
eq. (|7.1() now lead to (iterated) G* and S* exchanges. 

The matrix elements relevant to compute the width of Q* are reduced via 



X] i®*\'^l!mdma(3D^(i)aRf3\{KN)j^^^ 
E (0*l«S!L^™"/3^7(O«^8/3|(A'A^)/=l, 



i=± 



(G*+ T \d3af3D+aDs^\n T) . (8.4) 



The relative factor l/\/2 in comparison with eq. 1)5.6(1 is essential so that for G* we have to replace 
Xe Xe* and Ye Yq* in eq. (jHH by 



Xq* 



and Fe. = 2^1^ {e*+ T {dsapD+^Dspln T) • (8.5) 



In the flavor symmetric case we find 
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Thus, disregarding phase space effects (due to G* lying roughly lOOMeV above G"*" for Nc = 3), 
we obtain 



0* 



Xe 
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From eq. H8.6() we notice that in the large-A'^c' limit the factor | that appears in eq. 1)8. 3() is 
recovered. This is precisely what is required in order to obtain the BSA result in this limit. We 
have numerically verified that in the symmetry breaking case 

lim ^= lim ^ = 1 (8.8) 

for any value itik > wItt such that Tq* = ^Tq in that case as well, see also table Furthermore 
from eq. ()8.7|1 we also notice that for Nq = 3 the Q* width is even more suppressed in comparison 
to 0"*" (by roughly a factor of five in the symmetric case.) Finite symmetry breaking lowers the 
prediction for the 0* width even further, cf. figure [7| Thus we expect a very sharp B* resonance 
with a width of roughly lOMeV about lOOMeV above the pentaquark. This is opposite to 
the above criticized scenario based on matrix elements of the axial current to compute widths for 
hadronic decays, e.g. the authors of ref. [ij] obtained Fe./Fe = 3. ..6. 



IX. CONCLUSIONS 

We have thoroughly compared the bound state (BSA) and the collective coordinate approaches 
to chiral soliton models in flavor SU{3). For definiteness we have only considered the simplest 
version of the Skyrme model augmented by the Wess-Zumino and symmetry breaking terms. 
However, our analysis merely concerns the quantization of fluctuations about the classical soliton 
that generates baryons states. Therefore our qualitative results are valid for any chiral soliton 
model. 

Often the collective coordinate approach is identified with the rigid rotator approach (RRA) as 
only quantizing the spin-flavor orientation of the soliton. But, a sensible comparison with the BSA 
requires the consideration of harmonic oscillations in the collective coordinate approach as well. 
This generalizes to the rotation-vibration approach (RVA). Here we have studied that scenario 
exhaustively and in particular ensured that the introduction of such fluctuations does not double- 
count any degrees of freedom. Only in this way the model generates an interaction Hamiltonian 
that describes hadronic decays. In doing so, we have solved the long standing Yukawa problem 
in the kaon sector. Technically the derivation of this Hamiltonian is quite involved, however, the 
result is as simple as convincing: In the limit Nq oo, in which the BSA is undoubtedly correct, 
the RVA and BSA yield identical results for both the baryon spectrum as well as the kaon-nucleon 
S-matrix. This equivalence also holds when flavor symmetry breaking is included. In the first 
place this ensures that we have correctly introduced the collective coordinates. On top of that, 
this result is very encouraging as it clearly demonstrates that collective coordinate quantization 
is valid regardless of whether or not the respective modes are zero-modes. Though the large-A^c 
limit is helpful for testing the results of the RVA, we have also seen that taking only leading terms 
in the respective matrix elements is not trustworthy. 

There has been quite some confusion from the potential disagreements between the BSA and 
the RVA. There are essentially two reasons for that confusion: 

1) The BSA yields a slowly rising phase shift in the S = +1 channel that led to the misin- 
terpretation that no resonance was present. However, the scenario is that there is indeed a 
(broad) resonance which, unfortunately, is hidden by a repulsive background. 

2) The 0"*" width was incorrectly computed. Those computations did not even attempt to make 
contact with the BSA. 

Here we have resolved these puzzles: We observe a relatively sharp resonance in the RVA for 
Nc = 3 that evolves into the broad structure seen in the BSA as Nq increases. 
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Another major result of our calculation is that in the flavor symmetric case the interaction 
Hamiltonian contains only a single leading Nc structure of SU{3) matrix elements for the G"*" 
KN transition. Any additional structure only enters via flavor symmetry breaking terms in the 
Lagrangian. This is at odds with earlier approaches that essentially assumed any possible structure 
and fitted coefficients from a variety of hadronic decays under the assumption of SU (3) relations. 
We stress that this is not a valid procedure as already the embedding of the classical soliton 
breaks SU (3) and thus yields different structures for different hadronic transitions. In particular 
strangeness conserving and strangeness changing processes are not related to each other in chiral 
soliton model treatments. As for mx = rn-,^ only a single transition operator exists, there cannot 
occur any cancellation between different (leading Nc) structures to explain the probably small O"*" 
width. 

Though these general results are model independent, the numerical results that we obtain for 
the masses and the widths of pentaquarks are not. Moreover, the prediction for the latter may 
suffer from those sub leading I/Nq contributions that were not considered here. Our numerical 
estimates indicate that the is roughly 700MeV heavier than the nucleon with quite a small 
width of about 40MeV. We predict the 6* to be about lOOMeV heavier than the 0"*" with an 
even smaller width of the order of lOMeV. It should be kept in mind, that all quantitative results 
for masses and widths are model dependent and that their accurate predictions are very delicate. 
Despite of that, chiral soliton models do predict the existence of the low-lying pentaquarks Q and 
O* with strangeness, S = +1 and isospin, / = 0, 1. This is irrespective of details of the model 
and/or the adopted approach. 

Let us finally speculate about implications that indications for pentaquarks seen earlier in kaon- 
nucleon scattering 6^, 67 1 may have on chiral soliton models. At that time the structures observed 
in the 5 = +1 channel were assigned the quantum numbers Poi(1830), -Pi3(1810), -Do3(1790) and 
-Di5(2070). In chiral soliton models the latter can only be interpreted as a quadrupole excitation 
of a rotational ground state in the Pqi channel. The structure observed in that channel is actually 
higher in energy than both the P13 and Dq^ structures. It is thus suggestive that this Pqi structure 
should not be identified with the rotational ground state whose mass should thus be significantly 
lower than 1790MeV. This would be a further hint for a 0"*" pentaquark with Mq < 1700MeV. 
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Note added 

After this manuscript had been submitted, it has been commented on by T. Cohen |6^. We 
have rebutted his arguments in a reply 



APPENDIX A: CONSTRAINED EQUATIONS OF MOTION 

In this appendix we will derive the equation of motion 1)3.15(1 for the fluctuations in the RVA. 
We find it more convenient to perform that analysis in a general language for t^q, rather than for 
the P-wave projection, eq. H3.2() . In that language the kaonic zero mode wave-function carries 
two indices. They parameterize the rotation of kaon degrees of freedom (a = 4, . . . , 7) into any 
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direction (a = 1, . . . , 8) of SU (3), 



2U 



: sm— — Xifi, 



(Al) 



Only four of these wave-functions (6 = /? = 4, . . . , 7) are non-zero. They correspond to different 
choices of the zero-mode analogue of the complex isospinor in eq. (|3.2j) and are normalized such 
that J d^r Za{x)MK{r)zZ{x) = Sp-y where Mxir) is the metric function defined in eq. (|3.4|1 . The 
analogue of eq. (|3.7|) reads 



d-^r 2A(r) f8a(iz2{x)z^p{x) = ujq fs-fS 



(A2) 



We postpone the discussion of the rotation-vibration coupling terms involving the angular velocity 
to the end of this appendix and collect all other terms in the Lagrange function that are 0{N^) 
according to the rules, eq. 1)5.11) . The result of this tedious calculation is 



L 



-M-^ns + ^eKnl + ^ I d\ 



Mk (r) Va ix)r]c,ix)- r]a (x ) h^p (x ) r]p {x 




Ox^a I d^rz'j^ix) 



{m\ - ml) fsa^Dsa j d'^r z^{x)ri^{x) + O 



(A3) 



where dots denote derivatives with respect to time. Since we are mainly interested in relations 



between conjugate momenta, the explicit form of the hermitian potential h 



(that also contains 



the spatial derivative operators as indicated in eq. ()3.4j) ) is irrelevant. For the Skyrme model 
Lagrangian, eq. (|2.3|) it may be traced from the literature, e.g. refs. [2^. Is^]. The last two terms 
in eq. ()A3|) are new and describe the coupling between fluctuations and the collective coordinates. 
The very last term would be absent in a flavor symmetric world. In that case only a single SU{3) 
operator, that is linear in the kaonic angular velocity 0,^, couples the fluctuations to the collective 
coordinates. We will soon see that in the contribution proportional to in eq. HA3|) only the 
piece that stems from the Wess-Zumino term survives because of the constraints. 
The conjugate momenta of the kaon degrees of freedom are 



Ra 

n„(x) 



dL 
6L 



-QK^a + ^QK I d^rz'^ix) 



MK{r)r]f^{x) - —X{r)fsi3yr]^{x) 



(A4) 



Sr]^ix) 

Obviously, they are linearly dependent. 



^pMK{r)zi{x) + MK{r)^^{x) - —X{r)fs/s^?}^{x) . (A5) 



Ra 




d^r Zp{x) 



n/3(^) - -j=X{r)fs(3^r]^{x) 



(A6) 



The second piece stems from the Wess-Zumino term. It is unexpected because the conjugate 
momenta originate from parts in the Lagrangian that contain time derivatives of the chiral field, U. 
In turn the linear dependence between the conjugate momenta results from the connection between 
^'^n and ^!^^^/^} causing a linear dependence between i?„ and n„ only. However this argument 
only holds for local pieces of the action while the Wess-Zumino term is non-local. 
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We disentangle the collective and fluctuation pieces from the rj momenta 



X 



QK^f3MK{r)z^{x) and n«(x ) = Mx(r)7?„(x ; 



V3 



X{r)fsap?jp{x) (A7) 



and demand that the momenta conjugate to the collective coordinates do not contain any fluctu- 
ation parts, 



The linear relation, eq. I)A6|) then translates into the primary constraints 



n/3(x) - -j=X{r)fsi3yr]^{x) 



0. 



(A8) 



(A9) 



The corresponding secondary constraints require the fluctuations to be orthogonal to the zero-mode 



Xa = I z1{x) MK{r)?jf3{x) = . 



(AlO) 



These constraints are linear functionals of the fluctuations and their conjugate momenta. They 
satisfy the Poisson brackets {^a^Xp} = ^ab so that and r/^ are conjugate to each other in 
the constrained subspace [t^I- This is unaltered by the (unexpected) contribution to ^'q from the 
Wess-Zumino term. After projection onto the P-wave channel the constraints, eqs. ()A9|) and (|A10|1 
turn into constraints for the associated radial functions given in eq. ()3.12() . 

We need to find the equations of motion for Ha and rja and also to extract the Hamiltonian 
for the interaction between the collective and fluctuation degrees of freedom. We start by the 
Legendre transformation and add the constraints and Xa with Lagrange multipliers Qq and Pa 
respectively, 



H = -L- Rana+ / d^rUa{x)qaix) + aa^a + PaXc 



(All) 



We express this 0{N^) Hamiltonian in terms of the conjugate fields Ha and rja as well as the 
angular momenta Ra, 



H 



2] MK{r) 



na(x ) + — A(r)/8Q/3 Vpix) 



1 



+■ 



Ra 



d^rz^{x) 



^a{x) + --^X{r)fsa(3V(3{x) 
) fsalsDsa J d^r Z^{x) rj^{x) + PaXa ■ 



+ 2]^ '^^«(^)^a/3(^)^/3(^) 



(A12) 



The Hamiltonian obviously contains terms that are explicitly linear in the fluctuations or their 
conjugate momenta. These terms contribute to the Yukawa interaction, eq. (|7.1() . Additional linear 
terms arise from the Lagrange-multipliers which we compute in two steps from the equations of 
motion 



Uaix) 



6H 



1 



6Uaix) MKir) [ 
6H 



^aix) + --j=X{r)fsa/3r]f3{x) 



6r]a{x) 



^a{x) + —X{r)f8af3 V/six) 
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+ 



V3 

-(5pMK{r)zi{x) 



(A13) 



First, the constraints, eq. (|A9|) and HA10|) must be valid for all times. Therefore their derivatives 
with respect to time 



Xc 



(frz'p{x) 
d\z1{x)MK{r)?jp{x) 



n/3(^) - --^A(r)/8;37 











(A14) 



must also vanish. Second, we substitute the equations of motion (|A13|) into these relations and 
extract the Lagrange-multipliers, 



Rn 



Ola 



d^r \{r)z'^{x)hp^r]^{x) 



2 r 



fsapDsfs - [m]^ - ml) / d^r z^{x)r]p{x) , 



(A15) 



with ujQ and F defined in eqs. (|3.7|) and (|3.8j) . respectively. As the conjugate momenta, these 
Lagrange-multipliers obviously carry collective 



a 



(coll) 



7m "7! °^'"'^7©r"V9r 



(A16) 



as well as fluctuation pieces: = — Oq''"'' and Pa = Pa — Pa"^^^ ■ When substituted into the 
Hamiltonian, eq. ()A12|) the collective pieces provide the above mentioned additional contribution 
to the Yukawa interaction. In total the latter reads. 



2Rn 



hit 



4 
~3 



^fsp^ / d'rz'^ix) [2X{r)-u;oMKir)]v^ix 



QKDsahap j d^rz^ix) 



m\ 



3F 

89^ 



(A17) 



The formula, eq. ()7.ip displayed in the main part of the paper finally results from replacing the 
zero-mode wave-function according to eq. ()A1|) and the intrinsic fluctuations according to eq. (|5.4|) . 

We abstain from presenting further details on the homogeneous parts of the equations of motion 
for IIq, and rja- It completely suffices to note that by (i) projecting on the P-wave channel, 
(ii) omitting the explicitly inhomogeneous pieces proportional to Ra and D^a (flavor symmetry 
breaking) and (iii) substituting the fluctuation pieces for the Lagrange-multipliers: 
and Pa Pp, eqs. (|A13|) directly transform into eqs. (|3.13j) that we employed to compute the 
background phase shift in section III. In the main text the solutions of these homogeneous parts 
are denoted IIq and r]^, respectively. 

Finally we would like to comment on the treatment of the angular velocity Jls and its Nq 
scaling. To gain some insight, we compute the time derivative of Ra from [Ra, H] with H given in 
eq. 



dRg 



1 



1 



7^ I fiajjRiR/S - ( 7;; 1- ^^8 ) fsapRp 



Rs 



(A18) 
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While the first term, which is subject to ordering ambiguities, is suppressed by 0{l/\JNc)-, the 
second term, together with the constraint iJg = suggests that ilg should be counted as O (-/V^) 
and thus not be omitted as Nc — > oo. 

By isolating the contribution from the time derivative of the ansatz, eq. (|2.11|1 we find to 
linear order in the fluctuations, 

^^^^ = ±(^l=u, + ns^ ^3(^)y^^M [fsapKvp{x,t)] A3(t)t + O(^) . (A19) 

where we have used that [Ag, Uq{x)] = 0. Hence ilg may be absorbed in a suitable re-definition of 
the energy scale u; ^ uj' = to + ^-nd thus be omitted form the rotation-vibration coupling^^. 

as in eq. (|5.2j) . However, it may not be omitted from the "classical" level, eq (|4.1|) which does not 
contain the fluctuation energy to compensate for fig. 
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